A note on the radial solutions for the 
supercritical Henon equation* 

Vivina Barutello, Simone Secchi 

Dipartimento di Matematica ed Applicazioni, Universita di Milano-Bicocca. 
via R. Cozzi 53, 1-20125 Milano (Italy) 

Enrico Serra 

Dipartimento di Matematica, Universita di Milano. 
via C. Saldini 50, 1-20133 Milano (Italy) 

Mathematics subject classification: 35J60, 34B15. 

Keywords: Henon equation, Neumann problem, supercritical nonlinearity, 
radial solutions 



Abstract 

We prove the existence of a positive radial solution for the Henon equa- 
tion with arbitrary growth. The solution is found by means of a shooting 
method and turns out to be an increasing function of the radial vari- 
able. Some numerical experiments suggest the existence of many positive 
oscillating solutions. 



1 Introduction 

In 1982, W.-M. Ni wrote the first rigorous paper, [15], on an equation introduced 
ten years earlier by Henon in |14j as a model for mass distribution in spherically 
symmetric clusters of stars. This equation goes now under the name of Henon 
equation, and was originally coupled with Dirichlet boundary conditions: 

-Am = |a;|"uP in Bi, 
u>0, in Bi, (1) 

u — 0, on dBi, 

where Si = {x € | < l}, with iV > 3, a > and p > 1. 

The existence of solutions to ((!]) for p < = 2* — 1 is a standard exercise 
in critical point theory that can be solved by various simple approaches. 



*Partially supported by M.I.U.R., national project Metodi variazionali ed equazioni dif- 
ferenziali non lineari. 
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On the other hand, Ni's main resuh states that ([T]) has (at least) one solution 
provided p < 2* — l + ^^rj ^^^'^ thus enlarges considerably the range of solvability 
beyond the classical critical threshold p — 2* — 1. It is also simple to prove, by 
means of the Pohozaev identity, that for p> 2* — 1 + -^z^, problem ([T]) has no 
solution. 

Thus for the Dirichlet problem (U)) the picture is rather sharp: setting pa = 
2* + -^z^, problem ((!]) is solvable if and only if p < pa — 1. In the following we 
will call Pa the Henon critical exponent. 

The key observation in Ni's work is the fact that the presence of the weight 
a;|", which is radial and vanishes at a; = 0, allows one to gain compactness 
properties when one restricts the analysis to radial functions. Indeed, by means 
of the pointwise estimate ([15i) 

|^^(kl)| < C*-^^— f^ttj foi' almost every x G Bi, (2) 
\x\ — 

which holds true for any radial u £ Hq{Bi), one can easily prove that the 
embedding of rad(^i) dx) is compact precisely for p < — 1, 

and this is what one needs to prove existence of a solution to ([T|) . 

Quite recently, much attention has been devoted to symmetry-breaking is- 
sues, namely to the question of whether least-energy solutions of ([T]) are radially 
symmetric functions (for example for large a) . In their seminal papers |20l I21j , 
Smets et al. proved that this is indeed false for a sufficiently large. In the 
last few years various aspects of the Henon equation have been analyzed, and 
the resulting literature is nowadays rather rich (see for example [1] , [B] , [7j , [5] , 
[9], [IQ], [17], [H], [H] and references therein). All these papers concern the 
Dirichlet problem. 

The Henon equation just recently started to draw attention when coupled 
with Neumann boundary conditions. In this case the problem reads 

' -Au + u^\x\°'uP inSi, 

< u>0, in Bi, (3) 

^f^ = 0, ondBi, 

and has been studied in [13j , where the authors proved some symmetry-breaking 
results by connecting the question of symmetry of the ground states to the 
symmetry properties of extremal functions in some trace inequalities. In the 
paper [T3] of course it is assumed that p < 2* — 1, to use a variational approach 
in H^Bi). 

From the point of view of the mere existence of solutions for the Neumann 
problem Q , the situation presents both analogies and discrepancies with respect 
to the Dirichlet problem. Indeed, also in the Neumann case it is very easy to 
check that the problem admits at least one solution if p < 2* — 1, and it has 
been proved in [THj that Ni's result extends to ([3]): by using an version of 
inequality ([5]), it is simple to prove that ^ admits at least one (radial) solution 
for every p < pa ~ 1. 

If one wishes to complete the picture of the solvability for ([3]) as a function of 
p, like in the Dirichlet case, one has to face the fact that the Pohozaev identity 
gives no relevant information in presence of Neumann boundary conditions. To 
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our knowledge, it is not known whether the critical Hcnon exponent serves as a 
threshold between existence and nonexistence in ([S]). 

The purpose of this note is to fill this gap, by looking for radial solutions 
of ^ without any limitation on p. Of course without bounds on p we cannot 
make use of variational arguments, and we take instead an ODE viewpoint. Our 
main result is the following. 

Theorem 1. For every p > 1 and a > 0, problem ^ admits a strictly increas- 
ing radial solution. 

The preceding result shows that the uselessness of the Pohozaev argument 
for nonexistence of solutions is not a technical obstruction, but reflects a com- 
pletely different situation with respect to the Dirichlet problem. Though our 
arguments are rather simple, we point out that it is rather difficult to find in 
the literature existence results for elliptic equations without any growth condi- 
tion; some exceptions, for singularly perturbed elliptic problems can be found 
for example in [51 [3] . 

Our main result can be easily extended to the more general problem 

'-Au + w = in Si, 

< u > 0, in Bi, (4) 

On 

— =0, ondBi, 
y. av 

under suitable assumptions on <j) and /, but always with no growth restrictions 
on /. 

The last section of the paper contains some numerical experiments that 
suggest a rather surprising fact: there are choices of iV, a and p for which 
many (and possibly infinitely many) radial solutions of ([3]) seem to exist. We 
clearly state that this is only a numerical hint towards further research, since 
no rigorous proof has been written so far. By the way, we believe that for any 
choice of the parameters, problem ^ has exactly one radial, increasing solution, 
but we have to state this only as a 

Conjecture. For every > 3, a > andp > 1, the only solution of ^ which 
is positive, radial and increasing is that of Theorem [7J 

The matter of uniqueness for radially symmetric solutions of semilinear el- 
liptic equations on balls or annuli is a classical and often overwhelmingly diffi- 
cult issue. We refer to [T^] and the references therein for a short summary of 
known results. We have been unable to find uniqueness results in the literature, 
concerning nonautonomous equations like ([3]) with a monotonically increasing 
dependence on |a;|. 

2 An existence result for the Henon equation 

A radial solution for ^ must solve the ODE problem 

'-u"-^M' + M = r"uP, in (0,1), 
< u>0, in (0,1), (5) 

u'{0) = u'{l) 0. 
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We observe that u'{l) = corresponds to the Neumann boundary condition, 
while we require u'{0) = to obtain classical solutions. Since we do not impose 
any upper bound on p, variational techniques do not seem to be useful to prove 
any existence result. For this reason we use a shooting method, which consists 
in finding 7 > such that the solution u-y of the initial value problem 

'-u" - + u^r°'uP, in (0,1), 

< u>0, in (0, 1) (6) 

satisfies u'^{l) — 0. Equation can be written in the form 

{r^-^u')' ^r^-Hu-r"uP), (7) 

and it is quite natural to introduce the auxiliary function A{r) = r^^^u'{r). 
By the definition of A and ®, we deduce that A{0) = and that A increases 
strictly if and only if u{r) < c(r) where 

The curve c will play a crucial role in our discussion; its peculiarity is that it is 
asymptotic to the coordinate axes, i.e. 

lim c(r) — +00 and lim c(r) = 0^. (8) 

r — >0+ r — ^+00 

The next Lemma is well known, see for example 

Lemma 1. For every 7 > 0, problem ^ is uniquely solvable on [0, +00). Its 
solution Ury is continuously dijferentiable with respect to the initial value 7. 

We begin now a qualitative study of solutions to with the aim of proving 
the existence of (at least) an initial datum 7 > such that the corresponding 
Uj matches the Neumann boundary condition at r — 1. 

Take any 7 > 0, and consider the solution of The corresponding 
auxiliary function A.y{r) — r^~^u'^{r) is strictly positive from r = until 
intersects the curve c, see ([T]). After the first crossing, the derivative of A^ 
becomes negative. Therefore there exists a zero of A'^ , and the first intersection 
point 

r^ -.^ inf{r > | A'^{r) 0} inf{r > | Uy{r) = c(r)} 

between the graph of and that of c is well defined. By definition, u'^{r) > 
(at least) on the interval (0,r^). Let us call the first nontrivial stationary 
point of u^, namely 

inf{r > | u'^{r) ^ 0}. 
The next Lemma states in particular that R-y is well defined. 
Lemma 2. For every 7 > we have r-y < R^ < +00. 

Proof. Since the function A^ steadily increases on (0,r-y) and A{0) ~ 0, we 
deduce that A^(r^) > and then that u'^{r^) > 0. Hence r^ < R^. Let us show 
that R-y < +00. Assume, for the sake of contradiction, that R^ = +00. This 
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means that u'{r) > for all r > r^. Fix any Rq > (of course u'^{Ro) > 0) 
and choose S such that 



< <5 < 



This choice is possible since increases strictly and c decreases strictly, so that 
u^{Ro) > c(i?o)- Therefore 



r u: 
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(r) -1 > S > for every r > Rq. (9) 



We now integrate ©i on the interval [Rq, R], with R > Rq; we have 



= - «;(i?o) + (iV - 1) / -^dr 



_7V 

-Ro 



J?. ^J?, 
u-y{r)dr+ / r°'u^{r)dr 

> -<(i?o) - / %(r-) [1 - r"uP-i(r)] dr 
JRo 



and hence 



u^{r) [r"MT,(r)P~i - l] dr < ^^(Eo) < +oo. 



iRo 

for every R> Rq. However, thanks to we obtain 



R 

lim / u^{r)\r°'u^{r)P-^ -1] dr > lim {R - Ro)-fd = +oo, 

R~*+oc .1 fl— >+oo 



a contradiction that concludes the proof. □ 

Lemma 3. There results u'J^{R^) < 0. 

Proof. Indeed, from equation © we obtain 

-u'^iR^) = u^{R^) (1 - R^uP-\R^)) 

which is a positive quantity since u-y{R-y) > c{Rj). □ 

Lemma 4. There exists 6 > such that for every j < S there results R-y > 1. 

Proof. Lemma [1] implies that sup^gjp ^ |w-y(r)| < 1 for all 7 sufficiently small. 
Therefore lies below the curve c on [0,1], since c decreases and c(l) — 1 
independently of the parameters p and a. The claim follows from Lemma [2J □ 

Lemma 5. lim-y^+oo R^ — 0. 

Proof. For the sake of contradiction we suppose that there exist S > and a 
sequence (^k)k, Ik +00, such that R^^ > S for every k. Since u'^^ is strictly 
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positive on (0,(5/2), the function A^^ is strictly positive on {0,5/2]. Hence, for 
every k, 











"\''-^u,,{r){l-r'^u^-\r)) dr 



< fr) (1 - r"7rM dr, 



where the last inequality holds because u-y^, (r) > 7^ on [0, (5/2]. Since — > +cx), 

we can choose fcg such that, for every k > ko, there results 7^, ''^ < (5/2; for 
such values of k we can split the integral to obtain 



,-<p-l,/o 



dr 



0< jl >■''-■«-,.« 



dr 



Indeed M-y^ increases and the quantity (1 — r"7^ ^) is positive on the interval 
(0,7^7^''"^^^") and negative on the interval (7^7'"^"^''^", (5/2). Integrating we 
reach the contradiction 

1 7fe / 



< "(7.:^-^^/") (0 



N N + a \2 



< 



whenever jk > {^^^^^^^^ {d/2)^°'/^^^^\ namely for 7^ large enough. □ 

We are now ready to state our main result. 

Theorem 1. For every p > 1 and a > 0, problem ^ admits a strictly increas- 
ing radial solution. 

Proof. We have seen that the first critical point of is larger than 1 
for small values of 7, and smaller than 1 for large values of 7. Consider the 
map F: (0, +00) x (0, +00) R defined by F{'-f,r) = u'^{r). We have that 
F{j,R^) = 0, and Lemma [3] implies that 

d2F{j,R^)=u'!^{R^) <0. 

The Implicit Function Theorem shows that 7 1-^ R^ is a continuous and even 
differentiable function. Therefore, there exists 7 such that R^ = \. This means 
that is a radial solution of the Neumann problem □ 
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Remark 1. When p satisfies the condition 



p - I < 



v/(A^-2)2+4 - (TV - 2) 



a, 



(10) 



the Maximum Principle shows that every radial solution for the problem ^ 
can change its monotonicity at most once, since it can intersect the curve c at 
most twice. Indeed, computing 



c"{r)-- -c'{r)+c{r) = 

r 

(p-l)2r-2-a/(P-i) 



a{N - l)(p - 1) - a{a + p - 1) 



it is easy to see that since r £ [0, 1], condition ((TU]) implies that c is a subso- 
lution for the operator —A + /. Furthermore, a positive solution u for Q is a 
supersolution for the same operator. Therefore w = u — c is a supersolution that 
vanishes when u intersects c. From the Maximum Principle we deduce that w 
can vanish at most twice. 



3 A generalization 

Theorem [T] can be adapted to problem ^ by imposing some suitable assump- 
tions of the functions (j> and /. First of all, both functions <j) and / are defined 
(and continuous) on [0, +cxd), since we are looking for positive radial solutions; 
furthermore we require: 

(hi) (p is increasing, 0(0) = £ > and lim 0(r) = k G [/, +oo]; 

r — ^ + oo 

(h2) the function s ^ f{s)/s is strictly increasing; 
(h3) lim ^ = 1 (= +00 if ^ = 0); 

(h4) lim ^ = i (= if K = +oo). 

Under conditions (hl)-(h4) the equation u/f{u) = (j)(r) defines implicitly a 
continuous curve u — S^(r) which plays the same role as the curve c in the pre- 
vious section. Indeed if we call H{u) = u/f{u), then ^(r) = which 
decreases since decreases and increases. Furthermore ^ is asymptotic to 
the coordinate axes. Indeed from (h3) we get 

lim ^(r) = lim H^^ (0(?')) = limiJ"^(M) = +oo; 

similarly, from (h4), 

lim ^(r) = lim H'^ (0(r)) = lim H~^{u) 0. 

We can now proceed exactly as in Section [2] defining the auxiliary function A 
and the points r^, R^. Lemmas [51 and [S] can be proved also in this setting 
with minor changes. We then obtain the required generalization: 
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Theorem 2. Let (j), f : [0, +00) —^Mbe continuous functions. If assumptions 
(hl)-(h4) are satisfied then problem ^ admits a strictly increasing radial so- 
lution. 

Remark 2. Our assumptions are clearly satisfied by nonlinearities with arbi- 
trarily fast growth at infinity, like /(s) = exp(s) — 1, or f{s) = exp(7s') — 1 for 
7 > and q > 1. This latter case is particularly interesting, because it corre- 
sponds to Trudinger-Moser type problems without any restriction on q and 7. 
Though in this paper we work in dimension A'^ > 3, it is immediate to check 
that our results hold also for N = 2, which is the case of the Trudinger-Moser 
problem. 

Of course, (h3) requires £ — 0, i.e. (f>{0) — 0. Similarly, a nonlinearity that 
is superlinear at zero forces k to be infinite. Also, we notice that homogeneity 
of (j) plays no role, as long as (f) satisfies the above assumptions; for example 
0(r) — r" + r^ , with a, f3 > 0, is an admissible function. 

Remark 3. It is proved in [Ij that in dimension iV > 3, for any p > 1 and for 
large values of R, the problem 



does not have any positive radial solution whose derivative changes sign. We 
could not find any similar statement for a nonlinearity like \x\"uP. 

4 Some numerical results 

This section is devoted to the description of some numerical experiments. Such 
results are purely numerical and non rigorous; the purpose of the authors is, on 
one hand, to give some examples of the existence result proved in Section[21 On 
the other hand we want to point out some features of interest in the behavior 
of the solutions for the shooting problem ^ when 7 diverges to -|-oo. These 
numerical experiments seem to indicate that the structure of the set of radial 
solutions of problem ^ is still far to be understood, and deserves further study. 

4.1 The monotone solution 

The monotone solution for the Henon equation corresponds to a choice of the 
parameter 7 such that the first maximum point i?-y of the solution of the shooting 
equation coincides with 1. 




for |a;| < R 
for |a;| = R 
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In the table on the right we have col- 
lected some values of 7, depending on 
N, a and p, for which \R^ - 1| < 10"^. 
The starred values of p are the Hcnon 
critical exponents, that is p — Pa — I- 
In the first four rows of the table, we 
fix the dimension N and the exponent 
a and we choose three values of p: sub- 
critical, critical and supercritical. We 
observe that 7 seems to be a decreasing 
function of p. 

In the last row we have investigated the 
behavior of 7 as a becomes larger and 
larger while p has a supercritical value. 
The numerical results seem to show that 
7 continues to lie near 1. This is prob- 
ably due to the fast convergence of 
to as 7 —> +00, see Lemma [H 



4.2 Numerical evidence of oscillating solutions 

We have considered so far the existence of a radial solution with Neumann 
boundary conditions with the first stationary point at 1. Recalling Ni's results 
about the existence of oscillating radial solutions for some elliptic problem on 
M" (see [T^) we now investigate the existence on non-monotone solutions for 
our problem. We address two natural questions. Can the shooting solution, 
which is defined on the whole interval [0, -l-oo), have stationary points different 
from Rjl Can we choose a, p, and 7 such that one of such points coincides with 
1? Although we do not have any rigorous proof of these facts, some numerical 
experiments show that the answers to these questions strongly depend on the 
parameters a and p, so that it seems very unlikely to obtain a single general 
result. 

Fix for instance = 4 and a — 5. The critical Henon exponent Pa ^ I is 
in this case p = 8. For such values of the parameters, the monotone solution 
corresponds to 7 w 1.034. If we compute the solution on a larger interval 
we observe that it oscillates. Let us call i?" the n-th stationary point of this 
solution. Lemma [5] suggest that, as 7 +00, i?" decreases also for 71 > 1; for 
this reason we increase 7 in order to obtain i?" ~ 1 for some n > 1. To obtain 
i?^ w 1 we need to choose 7 « 155, while for w 1, a value of 7 ~ 2584 will do 
(see Figure [TJ . These cases point towards the possibility of the coexistence of 
multiple radial positive solutions. Rather surprisingly this interesting behavior 
gently disappears when p becomes supercritical. Indeed the oscillations become 
less and less sharp as p increases and disappear when p is greater than some p 
(in the described case p w 16). This behavior is illustrated in Figure [5] and it 
suggests a strong difference with the autonomous case studied by Ni in [TO] . 
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Figure 1: the three plots represent the shooting solution to problem ©, cor- 
responding to different values of 7, when A'^ = 4, a = 5 and p = 8 (critical). 
In the horizontal axis the radial variable r varies to the interval [0, 10]. In the 
figure above-left we choose 7 — 1.034 and we obtain the first stationary point 
R-y ~ 1; above-right, when 7 « 155, the second stationary point approximates 
1. If 7 w 2584 the third stationary point satisfies « 1. In the second and 
third picture, the first maximum point R-y is not printed in order to obtain a 
reasonable scaling, and we just plot the solution for r > 0.2. 
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Figure 2: in the picture we compare the plots of the shooting solutions when 
A'' = 4, a = 5 and 7 = 1.034, when p varies from the critical Hcnon exponent 
p = 8 to some supercritical values. In the horizontal axis the radial variable r 
varies in the interval [0, 10]. When p is critical the numerical solution oscillates 
sharply; as p increases the oscillations become weaker and weaker and they 
disappear when p > 16. 
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